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ABSTRACT 

In this paper, the unsteady MHD free convective flow through a porous medium bounded by an infinite vertical 
porous plate with heat sink is investigated. Approximate solution for the velocity field, temperature field, mean skin 
friction and mean rate of heat transfer are obtained. The effects of various emerging parameters on the velocity field, 
temperature field, mean skin friction and mean rate of heat transfer are discussed through graphs and tables in detail. 
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1. INTRODUCTION 

The problems in heat transfer have so many applications in Engineering Sciences such as the design of cooling 
systems for motors, generators and transformers. Chemical engineers are concerned with the evaporation, condensation, 
heating and cooling of fluids. An understanding of the laws of heat flow is important to structures and to the architect in the 
design of buildings. To estimate the cost, the feasibility and the size of the equipment necessary to transfer a specific 
amount of heat in a given time, detailed heat transfer analysis must be made. The dimensions of boilers, heaters, 
refrigerators and heat exchangers depend not only on the amount of heat to be transmitted but rather on the rate at which 
heat is to be transmitted under given conditions. It follows that in almost every branch of engineering, heat transfer 
problems are encountered which cannot be solved by thermodynamic reasoning alone but require an analysis based on the 
science of heat transfer. 

In the recent years, the problem of free convective unsteady flow through a porous medium bounded by an infinite 
vertical porous plate has attracted the attention of number of researcher because of its possible application in design of 
steam displacement process in an oil recovery and various geothermal systems. The free convective unsteady flow past an 
infinite porous plate with heat source was studied by Pop and Soundalgekhar [1]. Free convective flow of a electrically 
conducting fluids past a semi-infinite flat plate under influence of a magnetic field has been studied by Gupta [2] and Singh 
and owling [3]. MHD free convective flow with Hall current in a porous medium for electro solution was analyzed by 
Sattar and Alam [4]. Sahoo et al. [5] have studied MHD unsteady free convective past an infinite vertical plate with 
constant suction and heat sink. 

The study of unsteady hydro magnetic free convection flow of viscous incompressible and electrically conducting 
fluid past an infinite vertical porous plate in presence of constant suction and heat absorbing sinks has been discussed by 
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Sahoo et al. [5]. Approximate solutions have been derived using multi-parameter perturbation technique. It is observed the 
increase in magnetic field strength decrease with increase in magnetic field strength. The hydro magnetic mixed convection 
flow of an incompressible viscous electrically conducting fluid and mass transfer over a vertical porous plate with constant 
heat flux embedded in a porous medium was investigated by Makinle [6]. Using the Boussinesq and boundary-layer 
approximation, the fluid equations for momentum, energy balance and concentration governing the problem are 
formulated. It was found that for positive values of the buoyancy parameters, the skin friction increased with increasing 
values of both the Eckert number and the magnetic field intensity parameter and decreased with increasing values of both 
the Schmidt number and permeability parameter. 

In view of these, in this paper we studied the MHD free convective unsteady flow through a porous medium 
bounded by an infinite vertical porous plate with heat sink. Approximate solution for the velocity field, temperature field, 
mean skin friction and mean rate of heat transfer are obtained. The effects of various emerging parameters on the velocity 
field, temperature field, mean skin friction and mean rate of heat transfer are discussed through graphs and tables in detail. 

2. NOMENCLATURE 



8 
P 

G 



: Acceleration due to gravity 

: Coefficient of volume expansion 

: Conductivity of the medium 



Ec = 



V' 



c p (K-r) 



: Eckert number 



B n 



Kinematic viscosity 
Magnetic induction 



Fluid density 



Gr = 



Grashoff number 



M = 



p 



7 '2 



Magnetic Parameter 



Permeability of the porous medium 



Pr = 



_MC P 



S = 



K 

45V 

v o 



Prandtl number 



Sink Strength parameter 
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c P 


: Specific heat at constant pressure 


V' 


. OUVLLUll VCHJClLy 


T' 


• T f 24 ITT n f 24 r ' 1 1 1 Iff 3 
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K 


: Thermal conductivity 


t 


: Time 


(u,v',Q) 


: Velocity of the fluid 


V 


: Velocity of the fluid along the y Direction 



3. MATHEMATICAL FORMULATION 

We consider the flow of an incompressible viscous and electrically conducting fluid through porous medium past 
an infinite vertical porous plate in the presence of constant suction and heat absorbing sinks. Let the X - axis be taken 
along the plate and y - axis normal to the plate. The fluid subjected to a constant transverse magnetic field of strength B 0 . 

Neglecting the induced magnetic field and Boussinesq's approximation, the equations governing the flow can be written 
as: 

Equation of Continuity: 

W A 

= 0 (3.1) 



dy 

That is V' = 7^ (constant) (3.2) 
Equation of Momentum: 

du Tr ,du' n/ . / x d 2 u oBl , V , 

+ V -^ = gjB(T -7 ) + V — °-u --u (3.3) 

at ay ay p k 

Equation Energy: 

dt dy pC p dy C p \dy ) 

On disregarding the Joulean heat dissipation. The boundary conditions are given by 

u '=o, v' = -vt, r=r;+£(r;-r:)^as/^o 

u -> 0, T -> T as / -> oo (3.5) 
Here U is the velocity of the fluid in the upward direction, V' is the velocity fluid along the y - axis, Vq is the 
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suction velocity (constant), t' is the time, g is the acceleration due to gravity, (3 is the coefficient of volume expansion, 
V is the kinematic viscosity, p is the fluid density, (7 is the electrical conductivity, B 0 is constant transverse magnetic 
field, k is the permeability of the porous medium, S' is the heat sink strength, T' is the temperature of the fluid near the 
Plate, C p is the specific heat at constant pressure and K is the thermal conductivity. 



Introducing the following non - dimensional variables 



y , 0 t'V' 2 4VQ) U jU , r 

y = — —, t = — — , co = — — , u =—,v = — ,T 
v 4v v' 2 y; p 



Avo) 



pC„ 45V 
Pr = ^, 5 = 



K 



'2 



G = , Ec = 

V? 



v r 



'2 



CAT'-T) 



-,M = 



aBO 2 



v V' v 2 1 

V = — D = = — 

V' 2 Vq W 2 Da 



(3.6) 



Where Pr, G, S, Ec and M are respectively the prandtl number, Grashof number, Sink strength, Eckert number and 
Hartmann number With the help of equation (3.5) and equation (3.6), equation (3.2) and (3.3) become 



1 du du d 2 u 

G.T H - Mu - Du 

4 dt dy dy 



(3.7) 



PvdT 0 dT d 2 T Pr arr n J du^ 

And Pr — = — T + — ST + Pr Ec 

u dt dy dy 4 



And the modified boundary conditions are 
u = 0,T = \ + ee im at y=0 
And w— >0, T^O as y -> a 



(3.8) 



(3.9) 



To solve equation (3.7) and (3.8), we assume £ to be very small and the velocity and temperature in the 
neighborhood of the plate as Solution of the Problem 



u(y,t) = u 0 (y) + ee im u l (y) and T(y,t) = T 0 (y) + £e iax T x (y) 



(3.10) 



Putting (3.10) in equation (3.7) and (3.8), equating harmonic and non-harmonic terms and neglecting the 

2 

co-efficient of £ , we get 
d 2 u n du, 



+ — 2— (M +D)u 0 =-GT 0 



dy 2 dy 



(3.11) 



d 2 u, du 



dy 2 + dy ' 



M+D + — k =~GT 
4 ) 



(3.12) 



Equations Hamic, Non-Hamic, Equation 
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d 2 T dT Pr 

^f+Pr^--— (i&syL = -2Pr£c 
dy 2 dy 4 1 
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(3.13) 



(3.14) 



Using multi parameter perturbation technique and assuming Ec « 1 , we write, 
u o = u oo + u 01 , T 0 = Tqq + Ec T 01 

Vl^ — ^10 ^11' ^1 — "^10 ^11 



(3.15) 



(3.16) 



Using equation (3.15) and (3.16) in the equations (3.11) - (3.14) and equating the coefficients of Ec and Ec 
only, we get the following sets of differential equations for u m ,u 0l ,u l0 ,u ll and T Q0 T 0l ,T w ,T n , 



u w + Ko - (AT + £>)w 00 - -GT m 
u oi + Ki ~(M + D)u m - -GT m 



M+D + — 

v 4 y 



M 10 — ^10 



M +Da + — 

v 4 



Andr 0 ;+Prr o ;+^r 00 =o 



L oo 



r 0 ;+Prr 0 ;+^r 01 =-Pr( M ;) 2 



l£ + Prl£-*(ifl>-S)r io =0 



Tu+PrTu ~(ie)-S)T u =-2Pr( M >; 0 ) 



(3.17) 
(3.18) 

(3.19) 
(3.20) 
(3.21) 
(3.22) 
(3.23) 
(3.24) 



Solving these differential equations (3.17) - (3.24) with aid of the corresponding boundary conditions and then 
substituting the values in the relations (3.15) and (3.16), we obtain the mean velocity u 0 and mean temperature T 0 as well 

as Wj,7[ as 



EcPrG 2 X 



(a 2 -a 2 -M - D) 



- + EcXAe- a * y 1 + - , 

(a 2 2 -a 2 -M -£>) (a z 2 -a 2 -M-D) 



-<hy 
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{ Y>rEcG 3 D 4 e _ 2aiy 



EcVxG 3 D K 



EcPxG 3 D. 



(a 2 2 -a 2 -M -Df~ ' ~ 2 ~ ** ^ 2 < ~ 2 



(a 2 2 -a 2 -M -D) 2 {a 2 2 -a 2 -M -D) 



(3.25) 



^ + EcX A 



a\-a k -(M + D + — ) 



te-w—F- 



al~a 6 -(M + D + — ) 



{l + 2E c PvX 2 X 3 }e 



a 6 y 



+Ec2GPxX 2 \_D n e- (a * +a * )y - D l2 e Ha * +at)y - D l3 e Ha * +a ° )y + D l4 e- (a2+a * )y ~\ 



(3.26) 



T n =n+ 



Pr EcXG 2 



>e' aiy — 



PxEcG 2 



(al-a 2 -M -Df {a 2 2 -a 2 -M-D) 



\D x e 2aiy +D 2 e~ 2a2y - D 3 e~ {ch+a ' )y } (3.27) 



T x =e~ a6y +Ec{2PrX 2 X 3 e- a6y -2Pr X^D^^ -D % e <a ^ )y - D 9 e <Ch+Ch)y -ZV^ 2 ^]} 

(3.28) 



Where 



a 2= I(Pr + Vp?^),a 4= I(lW ^^ 



l + Jl + 4 



a„ = • 



M + D + 

V 4y 



. A = 



a A 



,A = 



a 0 



Pr 9 '" 2 
4«:-2Pr</. + ^ 2 



4a 4 -2Pra 2 + 



PrS ' 



A = 



2a 2 a 4 



A 



(a 2 + a 4 ) 2 - Pr(a 2 + a 4 ) + 



£) - ! D = 

PrS ' 4 (4o 4 2 -2fl 4 -M-D)' 5 (4a 2 - 2a 2 - M - D) 



D 6 = 



A 



[(a 2 +a 4 ) 2 -(a 2 +a 4 )-M-£>] 



A = 



^8 



Pr 

(a 4 + o 8 ) 2 -Pr(o 4 +a g ) (ico-S) 



A = . 



a 4 a 6 



(a 4 +a 6 ) 2 -Pr(a 4 +a 6 ) — ^-(i(D—S) 



Pr ... 
4 



, A = 



2 Pr 
(a 2 +^ 8 ) -Pr(^ 2 +^ 8 ) (ico-S) 



Pr 

(a 2 +^ 6 ) 2 -Pr(a 2 +a 6 )- — (ico-S) 



(a 4 + ^ 8 ) 2 -(a 4 +a 8 )-(M + 
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D 



12 



D a 



(a 4 +a 6 ) 2 -(a 4 +a 6 )-(M + D+"—) 



(a 2 + a„) - (a 2 + cl) - (M + D + — ) 



D, 



10 



(a 2 + a 6 ) 2 - (a 2 + a 6 ) - (M + D + 



id), 
T 



X =D X +D 2 -D 3 . 



X 



PrG 3 X 



PrG 3 



1 (a 2 -a 2 -M -Df (a 2 2 -a 2 -M -D) 2 



{D 4+ D 5 -D 6 }. 



X 2 = 



[a 2 2 - a 2 - M - D)\ a\ - a 6 - M - D - — 



,X 3 =D 1 -D S -D 9 +D U 



And X A = 



2GPrX 2 X 3 



2 r, 10). 

a 2 6 -a 6 -(M +D + -) 



2GPtX 2 [D n -D l2 -D l3 + D l4 ] 



Skin Friction 

The skin friction at the plate in dimensionless form is given by 



T = 



= u' 0 (0) + £e ia *u' l (0) 



(3.29) 



(3.30) 



d y)y=o 

Rate of Heat Transfer 

The rate of heat transfer at the plate is given by 

*w=(ln) =W)+«^(°) 
v d y) y =o 

4. RESULTS AND DISCUSSIONS 

Figure 1-4 show the profiles of mean velocity and transient velocity. Figure 1 shows the effects of Prandtl number 
Pr and Hartmann number on the mean velocity U 0 . It is found that, the mean velocity is greater for mercury 

(Pr = 0.025) than that of electrolytic solution (Pr = 1.0) . Also, it is observed that the mean velocity decreases with 
an increase in Hartmann number M . The effect of Darcy number Da on the mean velocity is shown in Figure 2 with 
fixing G = 5, E c =0.001, CO = 5, S =-0.05, Pr = 1 andM =1. It is observed that the mean velocity is 
increases with increasing Darcy number Da . 



www.tjprc.org 



editor@tjprc.org 



82 



L. Hari Krishna & R. Ramakoteswara Rao 



0.5 




0 1 2 3 4 5 

Figure 1: Effect of Pr and M on Mean Velocity for G=5, E c =0.001, co=5, S=-0.05 and Z)a=0.1 
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Figure 2: Effect of Da on Mean Velocity for G=5, E c =0.001, o=5, S=-0.05, Pr=l, M=l 
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Figure 3: Effect of Pr and M on Transient Velocity for G=5, £ c =0.001, o=5, S=-0.05 and Z>a=0.1 
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Figure 4: Effect of Da on Transient Velocity for G=5, £ c =0.001, co=5, S=-0.05, Pr=l, M=l 
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Figure 3 represents the effects of Prandtl number Pr and Hartman number M on the transient velocity U . 

It is found that, the transient velocity U decreases with increasing Pr . Also, it is observed that, the transient velocity 

decreases with an increase in M . The effect of Darcy number Da on transient velocity U is shown in Figure 4 by fixing 

G = 5, E c = 0.001, CO = 5, S = -0.05, Pr = 1 andM = 1 . It is observed that, with the increases in Da increases 

the transient velocity. 

Figure 5 depicts the effect of Prandtl number Pr on the mean temperature T 0 . It is found that, the mean 
temperature decreases with increasing Pr . The effect of Prandtl number Pr on the transient temperature T for 
G = 5, E c = 0.001, 60=5, S = -0.05, M = 1 and Da = 0.1 is shown in Figure 6. It is observed that, the 
transient temperature decrease with increasing Pr. 




o A 1 1 1 — i — 

0 1 2 3 4 5 

y 



Figure 5: Effect of Pr on Mean Temperature T 0 for G=5, £ c =0.001, io=5, S=-0.05, M=l and Da=0.1 




Figure 6: Effect of Pr on Transient Temperature T for G=5, £ c =0.001, co=5, S=-0.05, M=l and Da=0.1 

Figure 7 depicts the effect of Hartmann number M on the mean skin friction T w . It is observed that, with 
increase in magnetic field strength decreases the mean skin friction. Figure 8 represents the effect of sink strength S on 
mean skin friction T for G = 5, E = 0.001, CO — 5, M = 1, Pr = 1 and Da = 0.1. It is observed that, the 

W J C 7 ' ' 

mean skin friction decreases with increasing sink strength S . 
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Figure 7: Effect of M on Mean Skin Friction T for G=5, £ c =0.001, co=5, S=-0.1, Pr=l and Da=0.1 




Figure 8: Effect of S on Mean Skin Friction for G=5, £,=0.001, o=5, S=-0.1, Pr=l and Da =0.1 
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Figure 9: Effect of Da on Man Skin Friction for G=5, £ c =0.001, co=5, M=l, Pr=l and 5=-0.1 
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Figure 10: Effect of Pr on Mean Rate of Heat Transfer for G=5, £ c =0.001, o=5, 5=-0.1 and Da=0.1 
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The effect of Darcy number Da on the mean skin friction T w for 

G = 5, E c — 0.001, CO — 5, M = 1, Pr = 1 and S = -0.1 is shown in Figure 9. It is observed that, the mean 

skin friction T , decreases with on increase in Da . Figure 10 shows the effect of Prandtl number Pr on the mean rate of 

heat transfer q w for G = 5, E c = 0.001, CO = 5, S = — 0. 1 and Da = 0. 1 . If is found that the mean rate of heat 

transfer decreases with increasing Pr . 
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Figure 11: Effect of 5 on Mean Rate of Heat Transfer q w for G=5, £ c =0.001, o=5, Pr=l, M=l and Z)a=0.1 
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Figure 12: Effect of Da on Mean Rate of Heat Transfer q w for G=5, £ c =0.001, co=5, Pr=l, S=-0.1 and M=l 

Figure 11 depicts the effect of sink strength S on the mean rate of heat transfer q w 
for G = 5, E c = 0.001, CO = 5, Pr = 1 M = 1 and Da = 0.1 . It is observed that, with increase in sink strength 
decreases the mean rate of heat transfer. The effect of Darcy number Da on the mean rate of heat transfer q w for 
G = 5, E c = 0.001, CO = 5, Pr = 1 S = -0. 1 andM = 1 is shown in Figure 12. It is observed that the mean rate 
of heat transfer q w increases with increasing Darcy number Da 
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Table 1: Table for Mean Temperature T 0 





G = 5 3 E € = 0.00L 
£> = 5 = S=-0.05 = 
andZ><3 = 0.1 


G = 5 = E C = 0.001= *>=5 = 
S=-0.05 s Pr = l andA/ = l. 


G = 5^E E = 0.001= 
Da = Q_\__Pr = \ . 
co — 5 and M = 1 


y 


M= 5 


M=30 


= infinite 




£k = Q.i 


S=-QA 


S=-0_2 


0 


1 


1 


1 


1 


1 


1 


1 


0.5 


0.6029 


0.602S 


0.6052 


0.605 


0.6029 


0.6029 


0.6034 


1 


0.3635 


0.3634 


0.3657 


0.3 65 5 


0.3635 


0.3635 


0.3 64 


1.5 


0.2191 


0.219 


0J2212 


0.22 11 


0.2192 


0.2192 


0.2195 




0.1321 


0.132 


0.133S 


0.1337 


0.1321 


0.1321 


0.13 24 


2 _5 


0.0796 


0.0796 


0.0S09 


O.OSOS 


0.0796 


0.0796 


0.079S 


3 


0.04S 


0.04S 


0.04S9 


0.04SS 


0.04S 


0.04S 


0.04S1 


3.5 


0.02S9 


0.02S9 


0.0295 


0.0295 


0.02S9 


0.02S9 


0.029 


4 


0.0174 


0.0174 


0.017S 


0.D17S 


0.0174 


0.0174 


0.0175 


4.5 


0.0105 


0.0105 


0.0107 


0.0107 


0.0105 


0.0105 


0.0105 




0.0063 


0.0063 


0.0065 


0.0065 


0.0063 


0.0063 


0.0064 



Table 2: Table for Mean Temperature T 





G = 5._E C =0.001 = 
0 = 5,5" =-0.05, 
andDa = 0.1 


G = 5._E C =0.001 = ^=5 = 

5 = -0.05 = ^ = l andM = l. 


G = 5._E S =0.00L 
Da =0.1 = Pr = l = 
0)-S and M - 1 


y 


M= 5 


M=30 


£>& = infinite 




£5=0.1 


S= - 0.1 


S=-02 


0 


1 


1 


1 


1 


1 


1 


1 


0.5 


0.5634 


0.5601 


0.5639 


0.5637 


0.5-617 


0.5632 


0.555-S 


1 


0.3252 


0322 


03253 


03251 


03232 


0325 


0.3162 


1.5 


0.1923 


0.1901 


0.1925 


0.1924 


0.1906 


0.1922 


0.1S43 


2 


0.1159 


0.1147 


0.1162 


0.1161 


0.1146 


0.1159 


0.109S 


25 


0.070S 


0.0702 


0.0711 


0.071 


0.0699 


0.070S 


0.0663 


3 


0.0435 


0.0433 


0.0437 


0.0436 


0.0429 


0.0435 


0.0404 


3.5- 


0.026S 


0.0267 


0.026S 


0.026S 


0.0263 


0.026S 


0.0246 


4 


0.0164 


0.0164 


0.0164 


0.0164 


0.0161 


0.0164 


0.0149 


4.5- 


0.01 


0.01 


0.01 


0.01 


0.009S 


0.01 


0.009 




0.0061 


0.0061 


0.0061 


0.006 


0.005-9 


0.0061 


0.005-4 



Table 1 Show the effects of Hartman numbers M , Darcy number Da and sink strength S . It is observed that, 

with the increase in Hartmann number reduces the mean temperature T 0 slightly. The mean temperature T 0 increases with 

decreasing sink strength S . Also it is observed that, the mean temperature T 0 increases with increasing Darcy 

number Da . Table 2 show the effects of Hartmann number M , Darcy number Da and sink strength S on the transient 
temperature T . It is found that the transient temperature T decreases with an increase in M . Also it is observed that, the 
transient temperature increases with increasing Da . Also, it is found that with increase in sink strength decreases the 
transient temperature. 
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